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ABSTRACT 

If at cosmological distances, a small fraction of gamma-ray bursts should be multiply 
imaged by intervening galaxies or clusters, resulting in the appearance of two very 
similar bursts from the same location with a relative time delay of hours to a year. 
We show that microlensing by individual stars in the lensing galaxy can smear out 
the light curves of the multiply imaged bursts on millisecond time scales. Therefore, 
in deciding whether two bursts are similar enough to qualify as multiple images, one 
must look at time scales longer than a few tens of milliseconds, since shorter time 
scales are possibly rendered dissimilar by microlensing. 
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1 INTRODUCTION 

After more than two decades of observations, Gamma-Ray 
Bursts (GRBs) remain arguably the most puzzling phenom¬ 
ena in astronomy. In particular, their distance scale is yet 
to be established to everyone’s satisfaction, with 100 kpc 
(Galactic corona) and 2 = 1 (~ few Gpc) being the main 
opposing candidates (see Fishman 1995). There are two ma¬ 
jor pieces of evidence that point to the cosmological origin 
of the GRBs. First, evidence of redshift in faint gamma-ray 
bursts has emerged, both from time dilation in faint bursts 
relative to bright ones (Norris et al. 1994, 1995) and from the 
observation that the peak energies of the spectra are lower 
for faint bursts (Malozzi et al. 1995), each effect placing the 
faint bursts at a redshift of about 1. Both effects are only 
detectable in a statistical sense because the durations and 
hardnesses of gamma-ray bursts range widely at any given 
brightness, and subtleties of data analysis therefore play a 
role, leading others to not detect the time dilation (Mitro¬ 
fanov et al. 1996). Also, there could be intrinsic correlations 
between the different properties of a gamma-ray burst that 
mimic redshift effects, though they are distinguishable in 
principle from the real thing, and in fact a cosmological ori¬ 
gin agrees better with the detected time dilation (see Wijers 
& Paczyhski 1994). 

Second, evidence has been found of correlations between 
bright gamma-ray bursts and objects at their expected red¬ 
shift (if cosmological) of a few tenths. Recently, Rood & 
Struble (1996) and Marani et al. (1996) showed that the 
brighter of the GRBs are correlated with Abell (AGO) clus¬ 
ters. Larson, McLean, and Becklin (1996) found an excess 


of galaxies at K band near well-localized gamma-ray bursts. 
The significance of these is not yet very great, but may im¬ 
prove with time. 

It has been pointed out that strong gravitational lens¬ 
ing, resulting in a GRB being split into more than one 
macroimage, may provide the best means to ascertain the 
cosmological origin of GRBs (Paczyriski 1986b). If a multi¬ 
ply imaged lensing event is observed, then the source can 
be assumed to be at a cosmological distance from us. How 
do we know if two bursts originating from the same loca¬ 
tion, to within the errors, are multiply imaged and are not 
due to a repeater or are simply two unrelated events? If the 
observations are free of noise, and microlensing by stars is 
unimportant, then the multiple images (macroimages) due 
to an intervening galaxy or cluster should have identical 
spectra and light curves, to within a scale factor equal to 
the relative image magnifications. The presence of noise will 
degrade the light curves, thus macroimages can appear dis¬ 
similar. The effect of noise and faintness of the macroimages 
on their potential classification as ‘identical’ was discussed 
by Wambsganss (1993) and Nowak & Grossman (1994). Here 
we consider microlensing by individual stars in the lensing 
galaxy, and show that it can ‘smear’ image light curves on 
time scales of up to a few tens of milliseconds, making mul¬ 
tiple images of the same burst even less similar. 

Should one expect to see multiply imaged GRBs if they 
are at cosmological distances? Mao (1992) and Grossman & 
Nowak (1994) estimated that the waiting time for one lensed 
pair to show up in the BATSE catalogue can vary from one 
to well over ten years, depending on cosmology and other 
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assumptions, so it is by no means certain that a lensed pair 
will be found with BATSE. 

Microlensing by stars is known to be important for some 
multiply imaged QSOs, for example, Q2237+0305, the Ein¬ 
stein Cross (Irwin et al. 1989, Houde & Racine 1994), and 
probably H1413-I-117, the Clover Leaf (Arnould et al. 1993), 
PG1115-I-080 (Vanderriest et al. 1986) and possibly a few 
more QSOs. Therefore it is expected that macroimages of 
GRBs will also be affected by microlensing. For concrete¬ 
ness, we examine one particular case that can hypothetically 
be observed. This case roughly corresponds to the image pa¬ 
rameters of a well known lens: Q2237-I-0305 (Wambsganss 
et al. 1990). We show what would have been observed if the 
source here were a GRB instead of a QSO. 

Throughout this paper we assume a standard cosmology 
with n = 1, and 
Ho = 75kms“^ Mpc“^. 


2 EXAMPLES OF MICROLENSED 

GAMMA-RAY BURSTS 

2.1 Microlensing calculations 

The gravitational potential of the leasing galaxy splits a 
single background burst into multiple macroimages, whose 
arrival times can differ by anywhere between a few hours to 
over a year, and whose angular separation will generally be 
less than several arcseconds. Because the galaxies are made 
up of stars, the leasing mass of the galaxy is grainy, and 
therefore the macroimages are also grainy, i.e. they consist 
of microimages that have a spatial as well as a temporal dis¬ 
tribution. Given the poor angular resolution of gamma-ray 
telescopes, not only the microimages, but also the macroim¬ 
ages cannot be precisely localized. On the other hand, the 
temporal resolution of some of the brighter GRBs is very 
good, of order milliseconds, and so the effects of microlens¬ 
ing are important here: the observed macroimages will ap¬ 
pear as a sum of many staggered and scaled versions of the 
original burst. 

To quantify the above description, we first need to spec¬ 
ify a macrolensing model, and hence the properties of the 
macroimages. Given these, we can obtain a temporal prohle 
of a microlensed 5-function burst, and later convolve it with 
an observed GRB light curve to produce an example of a 
possible microlensed GRB profile. 

We start with a macrolensing model of a galaxy, given 
by a non-singular isothermal sphere, with a core radius (See 
eq. [8.37] and section 8.1.5 of Schneider et al. 1992 for de¬ 
tails) . Let us now generate a set of macrolensing parameters 
reminiscent of the actual observed QSO lens, Q2237-I-0305. 
In that case, five images are observed, four of similar magni- 
hcation, and a demagnified central image. In our case, since 
the lens model is assumed to be circularly symmetric (for 
simplicity), only three images will be produced, two images 
of comparable magnification, and a demagnihed third one. 
Let the galaxy core radius be, rc=0.36 kpc {6c ~ 0.5"), 
and its velocity dispersion, (j„ = 185 km/s. The source and 
lens redshifts are Zs = 1.7, and zi — 0.04, respectively. For 
a source impact parameter, da = 0.125", three images are 
formed. The two brighter macroimages have magnihcations 
of 9.86 and 9.31, a separation of 1.72", and a time delay of 
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Figure 1. Microlensing light curves of a single <5-function burst. 
The macrolensing parameters are k=0.524, 7 = 0.354, with a total 
mean magnification of A = 9.86. Horizontal axis is the time delay 
in milliseconds, and vertical axis is the brightness of individual 
microimages, such that the total unmagnified brightness of the 
source is unity. Notice the diversity of microlensing profiles, and 
that the total time width generally exceeds several milliseconds. 


0.92 days. At the locations of these macroimages, the op¬ 
tical depth and shear (k, 7 ) due to the galaxy are (0.524, 
0.354), and (0.799, 0.385), respectively. These parameters 
are required to generate the microimage distributions. 

Given a set of At and 7 for the first of the two macroim¬ 
ages, we can now numerically simulate the action of mi¬ 
crolensing on a 5-function burst of radiation from a small 
source. The patch of the leasing galaxy where the macroim¬ 
age is formed is represented by a two-dimensional random 
distribution of stars. The mass function of stars is taken 
from Scalo (1986); it has an average mass of M = O.386M0, 
and an upper and lower mass cutoffs of 63Mq and O.O87M0 
respectively. A ray-tracing method implemented using a hi¬ 
erarchical tree code is used to find the microimages of a small 
source. Once the locations of the microimages are found, the 
geometrical and gravitational parts of the time delay are 
calculated for each microimage separately. Figure shows 
three representative examples of the resulting time series. 
Note that the three light curves look very different from each 
other, which is expected since microlensing at even moder¬ 
ate optical depths is a highly non-linear process. A good 
way to parameterize a given light curve is to calculate its 
brightness-weighted rms spread At, 


At = 


EL 


where t = 


EL 

EZi A. 


N is the total number of microimages in the lightcurve; tFs 
and Ai’s are their individual arrival times and magnifica¬ 
tions, respectively. It is evident from Figure ^ that the val¬ 
ues of At can vary appreciably, even for a fixed set of {k, 7 ). 


© 1996 RAS, MNRAS 000, 000-000 





















Microlensed GRB light curves 3 



Figure 2. Normalized cumulative histograms of time delay 
spread of microimages, Af in milliseconds, for a set macroimages 
with re=0.799, 7 = 0.385 (solid line), and k=0.524, 7 = 0.354 
(dashed line). 

Figure ^ is a normalized cumulative histogram of At values 
for many independent macroimages having the same k, 7 
values: the solid line is for (0.799, 0.385), and the dashed 
one for (0.524, 0.354). The median At is 5-10 milliseconds, 
and is therefore relevant for some GRBs. 


2.2 Application to real GRB light curves 

The magnitude and detectability of microlensing distortions 
depend both on the temporal structure and signal strength 
of the original burst. To illustrate this, we convolve two 
of the three microlensing patterns of Figure with two 
rather different bursts. For this purpose we constructed 
two time series that after binning and addition of appro¬ 
priate amounts of noise passably resemble BATSE bursts 
IB 910711 and 1B920218B (the latter shortened by a factor 
4 in time). The former was one of the shortest ever seen 
and contained a sub-millisecond spike (Bhat et al. 1992); 
the latter is more typical. Figure H shows a convolution of 
microlensing time profiles with the two simulated GRB light 
curves. As can be seen, the short burst can be distorted be¬ 
yond recognition, so it would take some courage to advance 
the view that they are a lensed pair (the paucity of such 
short bursts would help, of course). The second case is more 
promising, for its high signal to noise ratio allows one to both 
recognize the overall similarity at most time scales and the 
differences at the shortest ones (especially the width of the 
last spike). From a visual inspection of the first BATSE cat¬ 
alogue (Fishman et al. 1994) it appears that some 20 bursts 
have visible structure at time scales of 10 ms or shorter, out 
of 40-50 in which the plot allows one to identify such struc¬ 
ture, so the fraction of bursts in which signs of microlensing 
could be detected if they were a member of a lensed pair 
may be as high as 50%. 


3 TIME WIDTH OF THE MICROIMAGE 

DISTRIBUTION 

It is interestiw to note that the temporal extent of microim¬ 
ages (Figure pi) is much larger than the time delay between 
two images created by a single isolated star. In the latter 
case, time delay is approximately %GM/c?, which is about 
0.04 ms for a solar mass star-lens if the images are of com¬ 
parable brightness. This difference between the time width 
of microimages in a galaxy and those due to an isolated 
star occurs because a collection of stars, as in a galaxy, act 
together to produce a wider At than a single isolated star. 

Before we proceed to derive a rough analytic expres¬ 
sion for the total time spread of a microlensed burst, we 
note that our analysis is related to that of Katz, Balbus 
& Paczyriski (1986). However, there is a difference between 
the two approaches that makes Katz et al. result not di¬ 
rectly applicable to our case. Just like we do below, Katz 
et al. look at the distribution of microimage magnifications. 
They pose a question, “What is the size of the region in the 
lens plane that catches a given percentage of the macroim¬ 
age flux?”. Below, we answer a somewhat different question, 
“At what distance from the source does the magnification 
of individual microimages become a certain fraction of the 
macroimage flux?”. 

Imagine a two-dimensional star field of microlenses of 
optical depth ac, and let all stars have the same mass of 
one mass unit. Consider a background source whose unper¬ 
turbed position coincides with the centre of the coordinate 
system. Roughly, there will be at least A-l-1 images formed: 
a primary image, and one image associated with every star- 
microlens in the vicinity of the source. In other words, the 
macroimage will consist of a ‘cloud’ of microimages sur¬ 
rounding the unperturbed source position. (For a pictorial 
representation of a microlensing ‘cloud’ see Paczyhski 1986a, 
and Wambsganss 1990). The size of the cloud is directly re¬ 
lated to the geometrical time delay of the microimages. Since 
the gravitational and geometrical parts of the time delay are 
of the same order, one simply needs to know the size of the 
cloud in order to estimate the total time delay of microim¬ 
ages. We will now calculate the radius of this cloud. 

Gonsider a microlens on the x-axis, a distance L from 
the origin, where L is large compared to the Einstein ring 
radius of a microlens. An image due to this microlens will be 
located at x = L-\-x', where x' is small. There are two contri¬ 
butions to the deflection angle in this case: the microlens it¬ 
self, and a sheet of microlenses interior to L. On average, the 
mass of this sheet is kL^ . The total deflection angle is then 
given by Oj, = kL?/ x-\- 1/x', and = 0. The magnification 
of the microimage, A~^ = 1 — 7 *, is due to shear only^ since 
by assumption there is no continuously distributed matter. 
The shear is given by 7 * = k — x' ^ ^ x' ^. Solving the 
lensing equation, 0 = x — kI? jx — 1 /x', for the location of 
the images, one finds that 1/x' ~ (1 —«:)L. Therefore, we ar¬ 
rive at the following relation between the source-microlens 
separation and the magnification of the corresponding mi¬ 
croimage, G ^ A~^{1 — k)~*‘. 

So far we have assumed external shear to be 0, there- 

* Note that this shear, 7 *, is different from the external shear, 7 . 
The latter is due to the overall potential of the galaxy, while the 
former is due to the matter distribution close to the macroimage. 
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Figure 3. A short (top row) and long (bottom row) synthetic gamma-ray burst affected by microlensing. The leftmost panel shows the 
ideal burst. The next two panels show how it is changed by folding the ideal GRB with the top microlensing pattern of Figure before 
and after binning and adding realistic noise. The next panels are the same for folding with the bottom microlensing pattern of Figure hi. 


fore the distribution of microimages is circularly symmet¬ 
ric with the diameter proportional to (1 — If exter¬ 

nal shear is introduced, the overall shape of the microimage 
cloud will get elongated, with the major and minor axes pro¬ 
portional to (1 — K — 7 )“^ and (1 — K + respectively^ 
Since it is the longer of the axes that determines the to¬ 
tal width of microimage distribution, the final expression is, 
4-1(1-k-7)-4. 

The total width of the arrival times of various microim¬ 
ages is roughly the time delay of the outermost microim¬ 
age of sufficiently non-small A. In our dimensionless units, 
a single isolated microlens will produce a time delay of 1 
between microimages of comparable brightness. By compar¬ 
ison, a sheet of microlenses will have a width of time de¬ 
lays of ~ — K — 7 )-^. For example, if images as 

faint as 0.01 times the source brightness are considered, i.e. 
A = 0.01, and k = 0.524, and 7 = 0.354 (as before) the total 
temporal extent of microimages is 671 dimensionless units, 
or, accounting for the M and (1 -I- zi) time dilation factor, 
about 10.8 milliseconds. 

The straight line in Figure ^ is the analytical derivation 
presented above. The points are the results of our numerical 
simulations using the ray-tracing code. The y-coordinate of 
each point, (At), is the average of At, as defined at the end 
of Section 2.1, over many realizations of the source position. 
A range of n, 7 parameters were used to test the analytical 
formula. The fit seems to be good considering the simplicity 
of our derivation. 


t This is valid if re < 1. If k > 1, the major and minor axes are 
proportional to (1 — re -|- 7 )“^ and (1 — re — 7 )“^. 



Figure 4. Analytically-derived (straight line, see Section g) and 
numerically calculated (points) relation between the rms spread of 
microlensing delays and the macrolensing quantity, (1 — re — 7 )“^. 
Empty circles and squares correspond to the (re, 7 ) cases repre¬ 
sented by solid and dashed histograms, respectively, in Figure ^ 
The straight line is, (At) = (1-1-2;) (SGM/c^) M A~^l'^ (1 — re — 
7 )—^, where 2 ; = 0.04 (reminiscent of the 2237-1-0305 quadruple 
QSO lens), A = 0.01 (see text), and M = 0.386Mq. 


4 GRB OBSERVABLES AND MICROLENSING 

In general, all macroimages of a multiply images GRB will be 
affected by microlensing to some degree, since the quantity 
(1 — re — 7 )-^ of macroimages is usnally greater than a few. 
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However what one really wants to know, when faced with 
a pair of candidate macrolensing events, is how much that 
particular pair has been affected. We will now show that 
there is insufficient information in the two gamma-ray bursts 
to estimate this independently. 

A macroimage is strongly affected by microlensing if 
the temporal extent of microimages is sufficiently large. The 
parameter that determines the temporal spread of the mi¬ 
croimages is (1 — K — 7 )”^. Therefore, one would want to 
know this quantity for each of the images. 

Given a pair of suspected macroimages, two parameters 
can be observed: the ratio of the total fluxes of the GRBs, 
which is roughly the magnification ratio of the two macroim¬ 
ages, A 1 /A 2 , and their relative arrival time, AT. The filled 
dot on the left panel of Figure ^ is an example of such a 
set of observables. The lines correspond to log(Ai/A 2 ) vs. 
dimensionless AT of the brighter of the two images pro¬ 
duced by a non-singular isothermal (solid lines), and a de 
Vaucouleurs (dashed lines) lens model. The two lines of each 
model have different central mass densities; since the iden¬ 
tity of the galaxy lens will not be known in such cases, its 
central mass density will be unknown to within a factor of 
2-3. The four dots on the right panel of Figure show the 
corresponding (1 — k — 7 )“^ values of the primary image, in 
these four models. It is seen that a set of A 1 /A 2 and AT 
does not constrain the (1 — k — 7 )“^ of the macroimage. 

To make things worse, there is uncertainty in the val¬ 
ues of A 1 /A 2 and AT. The ratio of magnifications can be 
affected by microlensing, realistically, by up to a factor of 
2, or 0.3 in the log. The relative time delay is uncertain 
because it is proportional to [Dos/DoiDig)rQ, where D’s 
are the angular diameter distances between the source, lens 
and observer, and ro is the length scale appropriate to the 
macrolens model: core radius for the isothermal model, and 
effective radius for the de Vaucouleurs model. For realistic 
ranges of GRB source redshifts, 0.8-1.2, and lens redshift, 
0.2-0.4, there is a 30-50% uncertainty in the cosmological 
parameter value; and a factor of at least 2 uncertainty in 
the core/effective radius. Thus there is a factor of at least 
4-5 uncertainty in the value of AT, and our lack of knowl¬ 
edge about the particulars of the lensing galaxy makes the 
determination of(l — k — 7 )“^ virtually impossible. 

Finally, (1 —k — 7 )“^ is related to the average of At’s; in 
any particular case, At can assume a range of values, even 
if(l — K — 7 )“^is known precisely (see Figures and |^. 


5 SUMMARY 

It has already been pointed out (Wambsganss 1993) that 
macrolensed copies of a single GRB need not look alike due 
to noise, and faintness. In this paper we discuss another rea¬ 
son: microlensing due to stars in the lensing galaxy. Most 
of the models of GRBs predict the ‘emitting region’ to be 
much smaller than the Einstein ring radii of even the small¬ 
est lenses. The former is of the order of 10^®cm, while the 
latter is roughly 10^'*cm for a 1 O“®M 0 compact object. The 
microlensing light curve of a short flash originating from a 
compact region is a series of spikes; the spikes are due to 
individual microimages, each with a different arrival time. 
The width of their time distribution is a function of the op¬ 
tical depth and shear of the galaxy at the location of the 



AT/(const D„/D„,DJ (l-zc-y)-* 


Figure 5. Relative macroimage magnifications, log(Ai/A 2 ), vs. 
relative time delay, AT, and {1 — k — 7 )“^ of the primary image. 
Four macrolens models are considered: two non-singular isother¬ 
mal spheres, and two de Vaucouleurs profiles, each with a dif¬ 
ferent, but realistic, central surface mass densities. Notice that a 
single set of log(Ai/A 2 ) and AT (large dot on the left panel) can 
correspond to four very different (1 — k. — 'y)~^ (four dots on the 
right panel). This example demonstrates that GRB observables 
are insufficient to determine the extent to which macroimages 
have been affected by microlensing. 

macroimages, and is 5-10 milliseconds for typical values of 
(«:, 7 ) of the lensing galaxy, and a typical stellar mass func¬ 
tion, with an average mass of a few tenths of Mq. Therefore 
GRB light curves are expected to be affected on these time 
scales. Similarly, if GRB spectra are varying on these time 
scales, they too will be affected. 
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